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Mean-field dynamo theory suggests that turbulent convection in a rotating layer of 
electrically-conducting fluid produces a significant a-effect, which is one of the key in- 
gredients in any mean-field dynamo model. Provided that this a-effect operates more 
efficiently than (turbulent) magnetic diffusion, such a system should be capable of sus- 
taining a large-scale dynamo. However, in the Boussinesq model that was considered 



by Cattaneo fc Hughes] ( 2006 ) the dynamo produced small-scale, intermittent magnetic 
fields with no significant large-scale component. In this paper, we consider the compress- 
ible analogue of the rotating convective layer that was considered by [Cattaneo fc Hughes] 



(2006). Varying the horizontal scale of the computational domain, we investigate the 
dependence of the dynamo upon the rotation rate. Our simulations indicate that these 
turbulent compressible fiows can drive a small-scale dynamo but, even when the layer is 
rotating very rapidly (with a mid-layer Taylor number of Ta = 10*), we find no evidence 
for the generation of a significant large-scale component of the magnetic field on a dynam- 



ical timescale. Like [Cattaneo fc Hughes (20061, we measure a negligible (time-averaged) 



a-effect when a uniform horizontal magnetic field is imposed across the computational 
domain. Although the total horizontal magnetic flux is a conserved quantity in these 
simulations, the (depth-dependent) horizontally-averaged magnetic field always exhibits 
strong fluctuations. If these fluctuations are artificially suppressed within the code, we 
measure a significant mean electromotive force that is comparable to that found in re- 
lated calculations in which the a-effect is measured using the test-field method, even 
though we observe no large-scale dynamo action. 



1. Introduction 

Turbulence can amplify magnetic fluctuations through the random stretching of field 
lines. Dynamo action is said to occur when this process is efficient enough to outweigh 
the effects of ohmic dissipation. In a hydromagnetic dynamo of this type, the magnetic 
energy usually grows at scales of a similar size to (or smaller than) the characteristic scale 
of the fluid motion. This is usually referred to as a small-scale (or fluctuating) dynamo. 
However, many astrophysical objects are characterised by magnetic fields that vary on 
much larger scales. Explaining the existence of these coherent large-scale fields remains 
one the greatest challenges in dynamo theory. Studies of large-scale dynamo processes are 
often based upon mean-field theory, which is a turbulent closure theory describing the 

f Email address for correspondence: b.favier@damtp.cam.ac.uk 



2 



B. Favier and P. J. Bushhy 



evolution of the large-seale quantities in terms of the statistical properties of the small- 



scale flow ( jMoffatt,,1978, ,Krause fc Radler 1980). In mean- field theory, the evolution 
equation for the large-scale field is derived from the induction equation by decomposing 
the magnetic field into mean and fluctuating parts. Following this procedure, assuming 
that there is no mean flow, it is possible to show that 



d{B) 
dt 



(1.1) 



where rj is the magnetic diffusivity, (B) is the mean magnetic field (defined using a 
suitable average (.)) and S = {u x b) is the electromotive force (emf) arising from the 
interactions between the fluctuating velocity field u and fluctuating magnetic field b. 
Although several recent studies have considered a more complex dependence of £ upon 



the mean field and its spatial derivative (see, for example, Brandenburg et al. 20086 



Hubbard & Brandenburg 2009), this mean emf is usually expressed in the following 
form, 



(u X b) 



a. 



(1.2) 



where Uij and f3ijk are coefficients that parametrise the small-scale effects in the mean- 
field model. In isotropic turbulence, the term involving (3ijk corresponds to an additional 
(turbulent) diffusive term in the mean-field equation. This tends to inhibit the growth 
of a mean magnetic field. The (Bj) term, on the other hand, plays a crucial role 
in driving dynamo action in mean-field models, via a process that is usually known as 
the Qf-effect. However, for this process to work, the underlying turbulence must lack 
mirror symmetry. This is certainly the case when the flow possesses helicity. Indeed, 
large-scale magnetic field generation has been observed in simulations of helically-forced 



homogeneous isotropic turbulence (Brandenburg 2001 Alexakis et al. 2005), although 
it should be noted that Brandenburg ( 2001[ ) found that it was necessary (arguably as a 
consequence of magnetic helicity constraints) to evolve the simulations for a reasonable 
fraction of an ohmic time-scale in order to generate significant large-scale fields. 

Natural dynamos are invariably found in rotating systems, where the Coriolis force will 
tend to drive helical motions. On the basis of mean-field dynamo theory, many naturally- 
driven flows should therefore be capable of sustaining a large-scale dynamo. This is 



certainly the case in near-onset, rapidly-rotating Boussinesq convection (Stellmach &: 
Hansen|[2004 ) , where the (kinetic) Reynolds number is small enough to ensure that the 



flow is rather laminar. However, the situation is less clear cut when the flow is turbulent, 
which is probably the case in most natural dynamos. For example, [Cattaneo fc Hughes] 
( 2006 ) found that no significant large-scale magnetic fields are generated by moderately 



rotating, turbulent Boussinesq convection, despite the presence of helicity. In fact, as 
in the corresponding non-rotating system ( Meneguzzi fc Pouquet||1989 Cattaneo]|1999 ), 
they found a strongly intermittent magnetic field distribution that was mostly confined 
to small-scale structures. One important conclusion from this study is that the presence 
of helicity does not necessarily lead to a large-scale dynamo. More recently, however, 



Kapyla et al. (20096) did find significant large-scale magnetic fields in a more complex 



model of dynamo action in rapidly rotating turbulent compressible convection, with a 
convectively-stable layer lying beneath the convectively-unstable region. 

There are various possible explanations for the presence of large-scale dynamo action 



in the model of Kapyla et al. (20096). The authors suggest (for reasons that we de- 



scribe below) that this result can be attributed to the fact that their system is more 



rotationally-dominated than that considered by Cattaneo fc Hughes (2006). This is cer- 
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tainly plausible, however the increased rotation rate is not the only possible explanation. 
The absence of stratification in the Boussinesq approximation may also be crucial. In 
an analytic study of an idealised compressible model, Steenbeck et al. ( 1966 1 derived an 
expression for the a-effect that depends upon the local density gradient. Moreover, if 
we return to Parker's original picture of an a-effect due to the action of many cyclonic 
eddies upon a magnetic field (Parker 1955), it is not unreasonable to assume that an 
asymmetry between the convective upfiows and downflows due to effects of stratifica- 
tion, which causes upwellings to expand and downwellings to contract, will lead to lower 
levels of cancellation between neighbouring flux loops. This would enhance the contri- 
bution of these small-scale motions towards the a-effect. However, there can be local 
asymmetries in the upfiows and downfiows even in Boussinesq convection (where, due 
to the symmetries inherent in the Boussinesq approximation, half-domain averages are 
used to calculate the a-effect), so this picture may be over-simplistic. Either way, further 
study is needed in order to clarify whether or not the effects of stratification enhance the 
a-effect in turbulent convection, thus promoting large-scale dynamo action. Another key 
difference between the simulations of Cattaneo & Hughes (2006) and those carried out 
by Kapyla et al. |20096} is th e presence of a stable layer in the latter set of calculations. 
Although Tobias et al. ( 2008 ) found that the inclusion of a stable layer did not promote 
large-scale dynamo action in their Boussinesq convective dynamo model, we also cannot 
discount the possibility it somehow promotes the generation of large-scale magnetic fields 
in the more rapidly-rotating compressible model that was considered by [Kapyla et al. 



(20096 ). Finally, it should also be noted that vertical magnetic field boundary conditions 



were adopted in most of the calculations that were described by Kapyla et al. (20096), 



whilst Cattaneo & Hughes ( 2006 ) always assumed perfectly conducting upper and lower 
boundaries. This may be an important point, although Kapyla et al. (20096) did find a 
a significant large-scale magnetic field in a single calculation with perfectly conducting 
boundaries. So the choice of boundary conditions does not appear to be crucial in this 
respect. 

It is clear that an unambiguous determination of the relevant mean-field coefficients 
would enhance significantly our understanding of the properties of large-scale dynamo 
action in rotating convection. Making certain simplifying assumptions, such as the (so 
called) first order smoothing approximation, it is possible to write down analytic expres- 
sions for the mean-field coefficients in terms of the statistics of the underlying turbulent 
flow. Unfortunately, the necessary simplifying assumptions are inappropriate for most 
natural dynamos and therefore cannot be relied upon in the case of these convectively- 
driven flows. Many studies have therefore focused upon the direct computation of the 



mean-field coefficients from numerical models of convection (see, e.g. Ossendrijver et al. 
2001 Giesecke et aL|20"05 ). Different methods can be used, depending on the coefficients 



desired. It is possible to measure a volume-averaged value for a (or a half-volume aver- 
aged value in the case of Boussinesq convection) by imposing a uniform magnetic field 
across the domain and then measuring the resulting mean electromotive force. If the 
imposed field is strong enough to be dynamically-significant, this will clearly inffuence 
the results of this process (see e.g. Hubbard et al. 2009), so the measurement of the 



kinematic value for a requires a weak imposed field. The measurement of the turbulent 
diffusivity is more complicated as it requires spatial variations in the large-scale imposed 
field. Therefore this quantity is usually measured using the test-field method, which was 
originally developed to compute the mean-field coefficients for convection in a spherical 
sheU ( [Schrinner erZ][2005l [20071 ). Since then, this method has been used to determine 
the mean-field coefficients for many different types of ffow, from homogeneous turbulence 
( Brandenburg et al]|2008a Mitra "al]|2009 ) , to compressible convection (Kapyla et al. 
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2009a). The idea of the test-field method is to investigate the response of the system 
to a prescribed set of linearly independent magnetic field vectors that may vary in both 
space and time (see, e.g. Schrinner et a/. [2007 Hubbard fc Brandenburg|2009 ). The mean 
electromotive force can then by measured for each test field, which in turn can be used 
to determine the mean-field coefficients. 

In rotating Boussinesq convection, the imposed-field method suggests that the a-effect 
is strongly fluctuating, even below the threshold for small-scale dynamo action where 
any magnetic fluctuations are produced solely by interactions between the flow and the 



mean field (Cattaneo & Hughes 20061. Furthermore, in the calculations of Cattaneo &: 



Hughes ( 2006 1 , the mean value of this fluctuating quantity was found to be negligible 
when compared to the rms velocity of the flow. As we have already discussed, these 
authors failed to flnd a large-scale dynamo in this case. However, when applying the 
test-fleld method to their compressible model, Kapyla et al. (2009al) find a much larger, 



well-defined a-effect, the peak magnitude of which is only weakly-dependent upon the 
rotation rate. Kapyla et al. (2009a) also find that the test-fleld method suggests that the 



turbulent diffusivity decreases with increasing rotation rate, which led them to expect 
that large-scale dynamo action is possible only in the rapidly-rotating regime in their 
simulations. Although, as we have previously discussed, the effects of compressibility (or 
even the underlying stable layer) could be playing an important role in promoting large- 
scale magnetic fleld generation in the calculations of Kapyla et al. ( 20096J , it is intriguing 
that the test-field method measures a substantial a-effect even at a comparable rotation 
rate to that adopted by Cattaneo & Hughes (2006). It is worth noting here that similar 
issues arise in the case of dynamo action generated by convection with rotation and 



shear. Such a system can produce large-scale magnetic fields (Hughes & Proctor 2009 



Kapyla et aL[2010a ) but different conclusions have been reached regarding the question 
of whether or this can be described as an a — dynam o. jHughes et al. (2011 ) make the 
point that Hughes & Proctor (2009) and Kapyla et al. (2010a) actually consider models 
that differ in several important respects. However differences in the methods that were 



adopted for measuring the a-effcct, in particular whether or not resetting (Hubbard & 
Brandenburg||2009 ) was used, also probably contributed to the lack of agreement in this 



particular context. 



The findings of Cattaneo & Hughes (2006) and Kapyla et al. (20096) raise certain 
questions regarding the relative merits of the imposed-field and test-field methods. Cer- 
tainly if (.) represents a volume average across the whole convective domain, or the upper 
half of the domain in Boussinesq convection, then the imposed-field method yields an 
unambiguous determination of the required component of the a-effect (at least in the 
case where the imposed fleld is weak enough that it does not perturb the flow). This 
quantity is obtained by dividing the appropriate component of S by the magnitude of 
the imposed mean field, which is constant if (as is usually the case) the boundary condi- 
tions ensure that the imposed flux is a conserved quantity. The situation becomes more 
complicated if the required component of a is to be regarded as a function of depth, 
and it is certainly possible that the a-effect may be more efficient at certain depths 
than at others (as suggested by the calculations of Kapyla et al. 2009&). In this case 
it is necessary to consider horizontally-averaged quantities. As pointed out by |Kapyla| 
et al. (2010a), a complication arises when using the imposed-field method to determine a 



depth-dependent value for a. Even if the imposed magnetic flux is a conserved quantity 
the horizontally-averaged magnetic field will normally be a function of depth and time. 



Although it has been argued (Childress & Soward 1972 Cattaneo & Hughes 2006 Hughes 
et aL|[20TT ) that the anisotropic nature of the horizontal averaging process implies that 



the evolution of the mean magnetic fleld is independent of the turbulent magnetic dif- 
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fusivity, the depth-dependence of the mean-field does nevertheless make it more difficult 
to measure the a-effect. In particular, there is no reason why the instantaneous value of 
the mean horizontal field should not pass through zero somewhere within the domain. So 
simply dividing the measured mean emf at a given depth by the instantaneous value of 
the mean field may not even produce a finite value for the required component of a let 
alone a meaningful one. This complication does not affect the test-field method. However, 
even taking into account the possible depth-dependence of the a-effect, further study is 
needed in order to explain why the different methods adopted by |Cattaneo &: Hughes] 
(2006) and Kapyla et al. (20096) yield such different measurements for the a-effect (even 
at comparable rotation rates). 

The aim of this paper is to determine the conditions (if any) under which compressible 
convection in a rotating layer can generate a large-scale magnetic field. Furthermore, we 
shall investigate the issue of whether or not it is possible to understand the properties of 
dynamo action in this system using the framework of mean-field dynamo theory, focusing 
upon some of the issues surrounding the measurement of the a-effect in numerical simu- 
lations. The governing equations, boundary conditions and parameters, together with the 
numerical methods that are used, are discussed in the next section. Considerations about 
the global hydrodynamic properties of the resulting convective flows are presented in sec- 
tion [3) In section |4] we discuss the dynamo properties of these flows. Section [5] is devoted 
to the measurement of the a-effect. Finally, in Section [6j we present our conclusions. 



2. Model and method 

2.1. Model and governing equations 



The model setup is identical to that described by Favier & Bushby (2012). We consider 
the evolution of a plane-parallel layer of compressible fluid, bounded above and below by 
two impenetrable, stress- free surfaces, a distance d apart. The upper and lower bound- 
aries are held at fixed temperatures, Tq and To 4- AT respectively. Taking AT > implies 
that this layer is heated from below. The geometry of this layer is defined by a Cartesian 
grid, with x and y corresponding to the horizontal coordinates. The z-axis points ver- 
tically downwards, parallel to the constant gravitational acceleration g = gz. The layer 
is rotating about the z-axis, with a constant angular velocity ft = flz. The horizontal 
size of the fluid domain is defined by the aspect ratio A so that the fluid occupies the 
domain < z < d and < a;, y < Xd. Various physical properties of the fluid, namely the 
specific heats Cp and Cy, the shear viscosity fj,, the thermal conductivity K, the magnetic 
permeability no and the magnetic diffusivity ry, are assumed to be constant. The system 
of governing equations is non-dimensionalised using the scalings that are described in 



Bushby et al. (2008). Lengths are scaled by the depth of the layer d. The temperature 
T and the density p are scaled by their initial values at the upper surface, Tq and po 
respectively. The velocity u is scaled by the isothermal sound speed \/R<,Tq at the top 
of the layer, where Ri, — Cp — c„ is the gas constant. We adopt the same scaling for the 
Alfven speed, which implies that the magnetic field B is scaled by ^/^qPqR^,Tq. Finally, 
we scale time by an acoustic time-scale d/^/R^TQ. 

Having non-dimensionalised the system, the governing equations for compressible mag- 
netohydrodynamics can be expressed in the following forrrjf] 



f Note that there is a typographical error in equation (2.2) of 
term V • [puu) should be replaced by pu • Vit. 



Favier & Bushby 



(20121, the 
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dp 
at 



-u-Vp-pV-u (2.3) 



= -u-Vu - KaTaJ'^z-Ku--VP+- (VxB) + 0(m + l)i + — V-r (2.4) 
ot p p p 

^ = -M . VT - (7 - 1) TV • w + ^V^T + ^fcll (arV2 + Co|V x (2.5) 

ot p p ^ 

dB 

~dt 

where the magnetic field satisfies W ■ B ~ 0. The pressure, P, is determined by the 
equation of state for a perfect gas, P = pT, whilst the components of the rate of strain 
tensor are defined by 



V X {ux B- kCoV X B) , (2.6) 



dui duj 2^ duk 
dxj dxi 3 dxk 



Several non-dimensional parameters appear in these governing equations. The parameter 
9 = at/To is the dimensionless temperature difference across the layer, 7 — Cp/c^ is the 
ratio of specific heat capacities, whilst m = gd/R^,AT — 1 corresponds to the polytropic 
index. The dimensionless thermal diffusivity is given by k = K / dpoCp{Ri,TQY^^ , a = 
pCp/K is the Prandtl number, whilst Co — V^pPo/K is the ratio of the magnetic to the 
thermal diffusivity at the top of the layer. Finally, Tap is the standard Taylor number 
Toq = 4,PqQ? d'^ I p? , evaluated at the upper boundary. 

Some boundary conditions are required in order to complete the specification of the 
model. In the horizontal directions, all variables are assumed to be periodic. The upper 
and lower boundaries are assumed to be impermeable and stress-free, which implies that 
Ux,z = Uy,z — Uz — Q eX z — Q and z = 1. The thermal boundary conditions at these 
surfaces correspond to fixing T = 1 at z = and T=l-|-0atz = l. In most of the 



calculations that are described in this paper, we follow Cattaneo & Hughes (2006) in set- 
ting Bz — Bx^z — By^z = at z = and z — \. This corresponds to the assumption that 
the upper and lower boundaries are perfect electrical conductors. However, as previously 



noted, Kapyla et al. (20096) usually adopted boundary conditions corresponding to a 
vertical magnetic field. In order to test the sensitivity of large-scale dynamo action to 
changes in the magnetic boundary conditions, we also perform two additional simulations 
(which are identical in every other respect to the perfectly-conducting boundary cases 
R3a and R4a, see Table [l] below) in which the magnetic field is constrained to be vertical 
at the upper and lower boundaries. This corresponds to setting B^ = By ~ Bz^z = at 
z = and z = 1. 

2.2. Model parameters and initial conditions 

These governing equations have a simple equilibrium solution, corresponding to a hydro- 
static, polytropic layer: 

T =l + 9z, p= (l + ez)"" , u^O, B ^0. (2.8) 

For all of the convective flows that are described in this paper, we adopt this polytropic 
state (plus a small thermal perturbation) as an initial condition. Once fully-developed 
hydrodynamic convection is established, a seed magnetic field (with zero net flux across 
the domain) is introduced into the flow. 

There are many non-dimensional parameters in this system. However, it is not possible 
to survey the whole of parameter space, so a number of these parameters are held constant 
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Run A Resolution 9 Ra r^" ^a k Co Re Rm Net helicity 

Rla 4 256^x120 3 3 x 10^ 332 0.0055 0.13 157 485 10"^ 

R2a 4 256^x120 3 5 x lO""' 20 10^ 0.0044 0.125 145 465 -0.15 
R2b 8 512^x120 3 5 x lO""' 20 10^ 0.0044 0.12 142 484 -0.15 

R3a 4 256^x120 3 1.5 x 10*^ 2.3 lO'' 0.0025 0.08 150 750 -0.33 

R3b 4 256^x120 3 1.5 x 10^ 2.3 lO'' 0.0025 0.12 150 500 -0.33 

R3c 8 512^x120 3 1.5 x 10^ 2.3 lO'' 0.0025 0.08 150 750 -0.33 

R3d 4 256^x120 0.2 1.5 x 10^ 2.6 lO'^ 0.0001 0.2 162 800 -0.04 

R4a 2 256^ x 120 3 4 x 10^ 0.9 10** 0.0015 0.06 150 1000 -0.38 

Table 1. Summary of the parameter values for the dynamo simulations. As described in the 
text, A is the aspect ratio, S is a measure of the thermal stratification, Ra is the Rayleigh number 
{Rac denotes the critical value for the onset of convection), Ta is the mid-layer Taylor number, 
K is the dimensionless thermal diffusivity, Co is the ratio of the magnetic to thermal diffusivities 
at the top of the layer, Re is the global Reynolds number and Rm is the corresponding magnetic 
Reynolds number. The net helicity is defined to be the time-average of (u • V x u) / (|ix ■ V x 



throughout this study. We choose a value for the polytropic index of to = 1, whilst the 
ratio of specific heats is given by 7 = 5/3 (which is appropriate for a monatomic gas). 
These parameter choices ensure that the polytropic state (as given by Equation ( |2.8[ )) 
is superadiabatically-stratified. We also fix the Prandtl number to be cr = 1. All other 
parameters in the system are varied (some more so than others). These parameter choices 
are summarised in Table [l] The main aim of this study is to address the effect of rotation 
upon large-scale dynamo action in a compressible convective layer. We therefore consider 
four different values of the mid-layer Taylor number (which is defined by Ta = (1 -f- 
e/2f '^Tao): Ta = 0, 10^ 10^ and 10^, which we refer to as cases Rl, R2, R3 and 
R4 respectively. The non-rotating Rla simulation is a reference case. The R2 cases 
correspond to a moderately-rotating regime that is similar to that considered by |Cattaneo| 



& Hughes ( 2006 1 . The objective of these particular simulations is to determine whether or 



not the absence of a large-scale dynamo in the calculations of Cattaneo & Hughes ( 2006 1 
can be attributed to the lack of stratification in their Boussinesq model. Finally, the 
rapidly rotating simulations R3 and R4 correspond to the regime considered by |Kap3da| 



et al. (2009& ), although it should be stressed again that there are a number of differences 



between the models (the possible effects of which will be discussed further in Section 4). 
Our aim here is to test the claim of Kapyla et al. (2009&) that the absence of a large- 
scale magnetic field in previous calculations of dynamo action in rotating Boussinesq 
convection (e.g. Cattaneo & Hughes 2006 Tobias et al. 20081 could be solely due to 
the fact that the rotation rate was not high enough. In all calculations, except case 
R3d, the thermal stratification is fixed to be 6* = 3, which corresponds to a moderately 
stratified layer (containing approximately 2.7 pressure scale heights). Case R3d is a quasi- 
Boussinesq calculation {9 = 0.2) of a similar type to those considered by Cattaneo & 



Hughes ( 2006 1 , albeit at a higher rotation rate. For each of the different rotation rates, we 
vary the aspect ratio, A. The largest value of A in each case is determined by the available 
numerical resolution. As the Taylor number increases, the scale of the horizontal motions 
decreases, which allows us to reduce A whilst maintaining a separation of scales between 
the convective motions and the width of the computational domain. 

The other two parameters that need to be specified are k and ^o- Rather than specifying 
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K directly, it is often more convenient to specify a mid-layer Rayleigh number 

Ra=im. + 1- m7) (1 + e/2f"'-' (!!i±Il^ , (2.9) 

which is inversely proportional to k^, and measures the destabilising effect of the su- 
peradiabatic temperature gradient relative to the stabilising effect of diffusive processes. 
Rotation tends to increase the critical Rayleigh number, Rac, for the onset of convection 
( Chandrasekhar|1961 ). Therefore adopting the same Rayleigh number in each case would 



produce convective motions that are less vigorous in the rapidly rotating cases. If the 
Rayleigh number is not to be kept constant, there are two obvious choices that could be 



made. Following the approach that is described by Stellmach & Hansen (2004), we could 
choose to keep the quantity (Ra — Rac)/Rac constant for each case. However, that would 
produce exceptionally high Rayleigh numbers in the rapidly rotating cases, which would 
require very high numerical resolution, so was not something that was feasible here. In- 
stead, we adjusted the Rayleigh number in each case so that the resulting hydrodynamic 
convective flows all had a similar global Reynolds number. In rapidly-rotating convec- 



tion, large-scale vortices can appear at high Reynolds numbers ( jChan|2007 Kapyla et al. 



2011). Clearly such flows are undesirable if we wish to find dynamo-generated fields on 
scales that are much larger than the characteristic scale of the underlying flow. Therefore 
we aim to maximise the Reynolds number subject to the constraint that we remain below 
the threshold for the production of these large-scale vortical flows. Deflning t/rms to be 
the root mean square velocity and Pmid to be the mean density of the mid-layer of the 
domain we therefore aim for a global Reynolds number (based upon the depth of the 
domain, which equals unity in these dimensionless units) of Re = PmidUrrns / {ncr) ~ 150. 
For completeness, the value of the quantity {Ra — Rac)/Rac is also shown in Table [T] 
for each of the simulations considered. It is clear that even though the Rayleigh number 
is increased with the Taylor number, we are closer to onset in the high Taylor number 
cases than in the non-rotating reference case. For we consider different values for each 
choice of Ta and Ra. Dynamo action is expected only in the low C,q regime, so we gener- 
ally attempt to minimise Co subject to resolution constraints (although we also consider 
the effects of adopting a higher value for Co in case R3b). Defining the global magnetic 
Reynolds number to be Rm = CA-ms/(KCo), this low Co regime is equivalent to the high 
Rm regime. 

2.3. Numerical method 
The given set of equations is solved using a modified version of the mixed pseudo- 



spectral/flnite difference code that was originally described by Matthews et al. (1995). 
Due to periodicity in the horizontal direction, horizontal derivatives are computed in 
Fourier space using fast Fourier transforms. In the vertical direction, a fourth-order fi- 
nite differences scheme is used, adopting an upwind stencil for the advective terms. The 
time-stepping is performed by an explicit third-order Adams-Bashforth technique, with 
a variable time-step. The aspect ratio is varied from A = 4 to A = 8. The resolution 
goes up to 512 grid-points in each horizontal direction and 120 grid-points in the vertical 
direction. A poloidal-toroidal decomposition is used for the magnetic field in order to 
ensure that the field remains solenoidal. 



3. General hydrodynamical considerations 

In this section, we briefiy discuss some of the hydrodynamic properties of the unmagne- 
tised convective flows that are adopted as initial conditions for the dynamo calculations. 
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Figure 1. The temperature distribution near the upper boundary (z = 0.1). Brighter regions 
correspond to warmer fluid whereas cooler areas of fluid are represented by darker regions. From 
left to right, the mid-layer Taylor number is increasing, from Ta = to Ta^ 10^. From top to 
bottom, the aspect ratio is increasing, from A = 4 to A = 8. For the sake of clarity, the aspect 
ratio of the case R4 has been increased from A = 2 to A = 4 by plotting four copies of the 
domain side by side (which explains the apparent four-fold periodicity in the temperature plot). 



Figure [T] shows the Taylor number dependence of the temperature distribution in a hor- 
izontal plane near the upper boundary (z = 0.1), for each of the moderately-stratified 
{9 — 3) cases that have been considered in this paper. Let us first consider the non- 
rotating A = 4 case (Rla). The motion consists of a few convective cells, and there is no 
clear scale separation between the typical scale of these cells and the horizontal extent 
of the computational domain. As the rotation rate increases, the convective cells become 
smaller. However, even at Ta = 10^ (where the motion is clearly influenced by the effects 
of rotation) it is necessary to increase the aspect ratio to A = 8 in order to obtain the 
required scale separation to study large-scale magnetic field generation. When the layer 
is rotating more rapidly, a A = 4 domain is probably sufficient to give adequate scale 
separation in the Ta — 10^ case, although we also investigated the dynamo properties of 
the flow in a A = 8 domain just to be certain that this is indeed the case. In the most 
rapidly-rotating turbulent simulation (at Ta — 10^), adequate scale separation can be 
achieved with a A = 2 domain. 

We deflne the brackets (.)^ to mean a spatial average over the horizontal coordinates 
at a given depth z, whereas the brackets (.) (without the z subscript) are defined to be an 
average over all spatial coordinates. The (time-averaged) depth-dependence of the root 
mean square velocity, t/rmsl^) = \/ is shown in figure^a). In all cases Urmsiz) is a 
slowly decreasing function of depth in most of the domain but increases with z near to the 
lower boundary. For a given value of z, the higher the value of Ta, the lower the value of 
Urmsiz)- In other words the convection is less vigorous in the rapidly-rotating cases. The 
horizontal integral length scale, lo{z) (as defined by Favier & Bushby 2012), is plotted 
in figure [2]jb). In the rotating cases, R2a, R3a and R4a, this quantity is only weakly 
dependent upon depth. In the non-rotating simulation, Rla, the integral scale decreases 
away from the upper and lower boundaries, taking a minimum value at approximately 
z = 0.75. As is clearly indicated by figure [l] lo{z) is smaller in the rotating cases. Note 
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Figure 2. The depth-dependence of various time-averaged quantities for hydrodynamic convec- 
tion: (a) The root mean square velocity Urnis{z), (b) the integral length scale lo{z) and (c) the 
relative kinetic helicity 'H{z). Similar results are obtained in the larger aspect ratio cases. 



that due to the stratification in the layer, lo{z) and U^^s{z) are both asymmetric about 
the mid-plane of the domain, particularly in the Rla and R2a simulations. 

These plots raise some interesting questions regarding the extent to which these sim- 
ulations are "equivalent". As already mentioned in section 2.2 we choose the Rayleigh 
number in these calculations so that the global Reynolds number is kept constant. How- 
ever, the global Reynolds number is defined in terms of the layer depth rather than a 
characteristic scale of the flow. If we were to define a local Reynolds number in terms of 
lo{z), as discussed in Favier & Bushby (2012), we would conclude that the rotating cases 
are less turbulent than their non-rotating counterpart. However, unless we are prepared 
to consider flows at lower Rayleigh numbers, it is not feasible with existing numerical 
resources to match (e.g.) the mid-layer local Reynolds number in all cases. This is why 
we chose to match the global Reynolds numbers in the present study, deferring further 
consideration of this problem to a later paper. It is worth noting here that [Kapyla et al.\ 
(20096 ) chose to keep the Rayleigh number constant as they increase the Taylor number, 
which means that even the global Reynolds number decreases with increasing rotation 
rate. This reduction in the intensity of the turbulence at higher Taylor numbers may 
be significant, although it is difficult to say whether or not this would be beneficial for 
large-scale dynamo action in rapidly-rotating convection. Of course, the global magnetic 
Reynolds number is also defined in terms of the layer depth. Were we to define a local 
value that was proportional to the integral scale, this would also decrease with increasing 
rotation rate. This reduction in the effective local value of the magnetic Reynolds number 
explains why we are able to adequately resolve dynamo calculations at higher (global) 
values of Rm in the rapidly rotating cases, without increasing the numerical resolution 
(see Table [T]). The relative merits of local and global definitions for these parameters are 
discussed in more detail in Favier & Bushby (2012). We choose to quote only the global 
values here for ease of comparison with previous studies. 

Mean-field theory suggests that the presence of kinetic helicity should enable the flow 
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Figure 3. (a) Magnetic energy versus time (in units of the turnover time). The magnetic 
energy is normahsed by the kinetic energy in the saturated phase, (b) T{z) for the same set of 
simulations. 



to generate large-scale magnetic fields. We define the mean relative kinetic helicity by 

(u - V X u) , , 

n{z) = — ^ . 3.10 

The effect of the Coriolis force is to give a negative correlation between the vertical 
vorticity and vertical velocity near the upper boundary, and the opposite near the lower 
boundary. The time-average of 'H{z), for the four different values of the Taylor number, 
is shown in figure [2jc). As expected, the mean helicity is almost zero in the non-rotating 
case. For the three rotating cases that are considered in this paper, the flow is highly 
helical, with a maximum value of |^| « 0.54 for Ta = 10^. In Boussinesq convection (see. 



for example, Cattaneo & Hughes 20061, this helicity profile is antisymmetric about the 
mid-plane of the layer. Due to the effects of stratification, this is not the case here. In 
fact, in all of the rotating cases, the relative helicity changes sign at z « 0.6. A possible 
measure of the net helicity in the system is given by the ratio (u-V x u) / (|m-V x u\). 
We show in Table [l] the time-averaged value of this quantity for each simulation. Unlike 
the Boussinesq case, there is clearly significant net helicity in the moderately-stratified 
rotating calculations. 



4. Dynamo simulations 

In this section we describe the dynamo properties of these convective flows. For each 
case, the procedure is the same: a weak seed magnetic field with zero mean is intro- 
duced into a fully-developed, statistically-steady simulation of hydrodynamic compress- 
ible convection. For most of our calculations (i.e. those which adopt perfectly-conducting 
boundary conditions for the magnetic field), the initial spatial structure of the magnetic 
seed field is given by B =(Acos(fci?/), Acos{kix), 0), where ki is a wave number which 
is comparable to (or smaller than) the wave number at which the kinetic energy peaks. 
The seed field is different for the simulation with vertical boundary conditions. In this 
case, we set B =(0, 0, A cos{kix) cos{kiy)) . For the seed field, A is the amplitude of the 
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initial perturbation, chosen so that the typical ratio between the initial magnetic energy 
and the kinetic energy in the hydrodynamic state is about 10~^. This initial magnetic 
field is then evolved in time until the dynamo saturates in the nonlinear regime. 

The time evolution of the magnetic energy for various different simulations is shown in 
figure |3][a). In each case, the magnetic energy is normalised by the kinetic energy in the 
saturated state, whilst the time is normalised by the mean turnover time, l/C/mis (which 
is obviously different in each case). For Ta — and Ta — 10^, the magnetic energy 
saturates at a level that is somewhere between 4 and 9% of the kinetic energy, with the 
rotating simulation saturating at a slightly higher level than the non-rotating case (for 



more details see Favier & Bushby 2012). For Ta = 10^, the magnetic energy saturates 
at approximately 20% of the kinetic energy for case R3a and slightly below 10% for case 
R3b. Given that the magnetic Reynolds number is larger in case R3a than it is in case 
R3b, it is unsurprising that the dynamo saturates at a lower level in the latter case. In the 
most rapidly-rotating case (R4a), the magnetic energy saturates at approximately 30% 
of the kinetic energy. It is also apparent from figure [3j^ a) that there are fewer fluctuations 
in the magnetic energy curves in the more rapidly-rotating calculations. This is a due to 
the fact that there are a larger number of convective cells in the rotating simulations, 
which improves the quality of the spatial averages. 

Like most previous compressible calculations, all of these simulations have been evolved 
for a period of time that is substantially less than the ohmic dissipation timescale (based 
upon the depth of the layer). Ideally, longer integrations would have been carried out, 
but this restricted period of evolution is entirely due to computational constraints: Large 
spatial resolution is needed in order to resolve the magnetic structures at high Rm, whilst 
the time-step stability constraints in compressible calculations are much more restrictive 
than they are in comparable Boussinesq (or anelastic) calculations. Having said that, 
all of these dynamos have saturated in the nonlinear regime. Furthermore, in their most 



rapidly-rotating case, Kapyla et al. (20096) observe a large-scale field after approximately 
70 turnover times (note that our definition for the turnover time gives a value that is a 
factor of 2tt smaller than the definition adopted by Kapyla et aL||20096 ) which suggests 



that we have evolved these calculations for a long enough period of time to identify 
large-scale dynamo action if it is present. 

4.1. Energy in the mean magnetic field 

Before considering the properties of these dynamo-generated magnetic fields in more 
detail, we first address the issue of whether or not these dynamos have produced a 



significant mean field. Following Cattaneo & Hughes ( 2006 ) , we define the quantity 



r(^) = ^ , (4.11) 

which is the ratio of the energy in the mean field to the total magnetic energy at a 
given depth. Figure [3|^b) shows the time-average of this quantity versus depth for several 
different cases. The first point to make is that, in all cases, F is very small at all depths. 
Therefore, for all the cases under consideration here, the energy in the mean magnetic 
field is very small compared to the total magnetic energy. At least in that sense, none 
of our convective flows appear to be generating significant large-scale fields. It is also 
worth noting that the higher aspect ratio cases systematically exhibit smaller values of 
r, which indicates that the computations are not in a regime in which the results are 
completely independent of the box size. The fact that F decreases when the box size 
increases suggests that small aspect ratio simulations tend to overestimate the ability of 
rotating convection to generate a large-scale dynamo. Perhaps surprisingly, the largest 
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Figure 4. Grey-scale plots of the horizontal component of the magnetic field, B^, a.t z = 0.1 
for (a) Ta — 10^ (run R2b) and (b) Ta — 10^ (run R3c). Contours are evenly spaced between 
Bx ~ —0.4 and Bx ~ 0.4. Light and dark tones correspond to opposite polarities. 



measured value for F corresponds to the non-rotating simulation (Rla). This simulation 
is also the simulation with the worst scale-separation (see figure [T]). Furthermore, at 
this aspect ratio (i.e. for cases Rla, R2a and R3a), the proportion of the total magnetic 
energy in the mean-field tends to decrease as the rotation rate increases. A similar trend 
is observed in the larger aspect ratio (A — 8) domain in cases R2b and R3c. Comparing 
cases R3a and R3b, which only differ by the value of the parameter Co, we see that they 
have a similar profile for F. This suggests that the relative strength of the mean field 
in these dynamo calculations is not strongly dependent upon the magnetic Reynolds 
number (at least for computationally accessible values of this parameter). Although the 
mean field is still negligible, it is interesting to note that the F is larger at all depths in 
the R4a case than it is in any of the R3 cases. Again, this effect could be attributed to the 
lower aspect ratio that has been used for this simulation. Although not shown here, our 
quasi-Boussinesq rapidly-rotating dynamo simulation (R3d) also exhibits a small value 
for F(z), with a typical value of F « 5 x 10"''. Therefore these findings do not appear to 
be sensitive to the level of stratification within the domain, a comment that is further 



reinforced by the fact that Cattaneo & Hughes (2006) also found small values for F(z) in 
their Boussinesq model. Finally, it is worth noting that the simulations that we carried 
out with a vertical magnetic field boundary condition (which were otherwise identical 
to cases R3a and R4a) also produced very small values for F(z), with typical values of 
F w 10"^ and F « 3 x 10"^ respectively. 

A visual inspection of the magnetic field distribution confirms, at least in qualitative 
terms, the absence of a significant large-scale field. Figure |4] shows the horizontal dis- 
tribution of Bx in the z — 0.1 plane, for both of the A = 8 calculations, approximately 
20 convective turnover times after the magnetic energy has saturated in the nonlinear 
regime. In both cases, the magnetic field distribution is dominated by small-scale features, 
with no obvious large-scale structures. We have also analysed the time-dependence of the 
horizontally-averaged horizontal magnetic field components (< >z and < By >z) as 



a function of depth, in a similar fashion to Kapyla et al. (200961 (see for example their 
figures 4 and 7). Figure [s] shows the time-evolution of these quantities (normalised by the 
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Figure 5. Horizontally averaged magnetic field components (a) < >z and (b) < By >z, 
as a function of depth, normalised by the equipartition magnetic field i?eq, for simulation R4a 
with vertical magneti c field at the boundaries. This can be compared with figures 4 and 7 from 



Kapyla et al. (200961 



1 /2 

volume-averaged equipartition field i?eq = {p\u\^) ) in the rapidly-rotating R4a case 
with vertical magnetic field boundary conditions, which is the case that is most compa- 
rable to the results presented by Kapyla et al. ( [20096 1. Unlike Kapyla et al. (200961, we 
observe no coherent horizontally-averaged magnetic fields. Similar results are obtained 
at lower rotation rates or with perfectly boundary conditions. 

4.2. Energy spectra for the turbulent calculations 

As discussed in the previous subsection, the quantity T{z) can be used to determine 
whether or not a dynamo is producing a significant mean horizontal magnetic field. 
However, it provides us with no information about any of the other large-scale magnetic 
modes. To gain further insight into the properties of these hydromagnetic dynamos, it 
is useful to consider their magnetic and kinetic energy spectra. For each horizontal wave 
number, fcjj, we define the horizontal magnetic energy spectrum in the following way: 



EhiikH) = ^ B{ka;,ky,z) ■ B*{kx,ky,z) 



(4.12) 



where B{kx, ky, z) is the two-dimensional Fourier transform of B{x, y, z) (with complex 
conjugate B* (kx, ky, z)) and the second of the summations is over all k^ and ky such 
that k'^ + ky = fc|f . Similarly, the horizontal kinetic energy spectrum is defined to be 



EK{kH)^j'^'^u{kx,ky,z) ■pu*{kx,ky,z) + u*{kx,ky,z) ■ pu{kx,ky,z) 



(4.13) 



Although we only discuss vertically-averaged spectra (as defined above), similar results 
are obtained when these spectra are computed at specific depths. 

We focus initially upon the spectra for the simulations that were carried out in the 
A = 8 domain. Figure |6](a) shows the kinetic and magnetic energy spectra, for both the 
kinematic phase and the saturated nonlinear phase, for the moderately-rotating Ta = 10^ 
case R2b, whilst figure [6]jb) shows the corresponding spectra for case R3c {Ta — 10^). 
During the nonlinear phase, both the kinetic and the magnetic energy spectra are time- 
averaged. In the more rapidly-rotating case, R3c, the kinetic energy spectrum peaks at 
kn = 15, whereas the corresponding kinetic energy spectrum in the R2b case peaks 
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Figure 6. Horizontal energy spectra for (a) case R2b (Ta = 10^ A = 8), (b) case R3c {Ta = 10 , 
A = 8), (c) case R4a {Ta = 10*, A = 2) and (d) R3a {Ta = lO'', A = 4). Due to the variations 
in the aspect ratio, the range of available horizontal wave numbers is different in each case. 
The magnetic energy spectra are plotted at different times during the kinematic growth phase 
(dotted thin lines) and time-averaged during the saturated nonlinear phase (dotted thick line). 
The kinetic energy spectra are averaged over time during both the kinematic and the nonlinear 
phases. The dash-dotted line on figures (a) and (b) is an indication of the power-law k^j^ at small 
wave numbers. 



at kn = 6. This simply reflects the fact that the horizontal scales of motion tend to 
be smaller at higher rotation rates. In both cases, these peaks in the kinetic energy 
spectra are at sufficiently large wave numbers that it is meaningful to talk about a 
separation in scales between the size of the convective cells and the horizontal extent 
of the computational domain. The kinetic energy spectra during the kinematic and the 
nonlinear phases are nearly indistinguishable on this logarithmic scale. However, for Ta = 
10^ (R.2b), there is some indication that the kinetic energy spectrum in the nonlinear 
phase is slightly reduced for modes in the range kn < 2. In the Ta = 10'' (R3c) case, 
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a similar reduction in the kinetic energy spectrum is observed for certain wave numbers 
in the range ku < 10. Note however that this reduction of the kinetic energy at low 
wave numbers during the saturated phase only corresponds to a very small fraction of 
the total kinetic energy, so no significant qualitative changes are observed in the flow as 
the dynamo saturates. 

In the A = 8 calculations, the spikes in the magnetic energy spectra (at = 6 
for the R2b calculation and k^ — 20 for R3c) during the first part of the kinematic 
phase correspond to the initial seed magnetic field. However, these magnetic fluctuations 
rapidly spread to all available scales. After a brief transient phase, the magnetic energy 
spectra appear to evolve in a self-similar way during the kinematic phase. Indeed, other 
than having a lower amplitude, the magnetic energy spectra during the latter part of 
the kinematic phase are almost indistinguishable from the final magnetic energy spectra 
in the nonlinear saturated phase. Ignoring the seed field, the magnetic energy spectra 
always peak at a higher wave number than the corresponding kinematic energy spectra. 
During the saturated phase, the wave number of maximum magnetic energy is k^ = 10 
for R2b [Ta = 10^) and fc// = 22 for R3c (To = lO'^). This behaviour is typical of 
calculations of small-scale dynamo action when the magnetic Prandtl number, Pm (i.e. 
the ratio of the global magnetic Reynolds number to the global kinetic Reynolds number), 
is greater than unity. This is the case in all of these dynamo calculations. The slopes 
that are indicated on figures [6](a) and [6](b) correspond to a scaling fc^ at low wave 
numbers. Given the anisotropy and inhomogeneity of the fiow (and the fact that only two- 
dimensional Fourier transforms are used here) we do not expect to recover the Kazantsev 
scaling, corresponding to /3 = 3/2. We find /3 = 1.1 ± 0.2 for case R2b {Ta = 10^) and 
/3 = 1 ± 0.1 for case R3c (Ta — Note that these values increase slightly if we 

consider only the middle of the convective layer (i.e. if the spatial average is restricted 
to 0.2 < z < 0.8), where the magnetic field is more isotropic. Our values are slightly 
larger than the ones predicted by Pietarila-Graham et al. (2010) and smaller than the 
ones predicted by Moll et al. (2011). However, this discrepancy is unsurprising. Their 
non-rotating calculations include many physical processes relevant to solar convection, 
and they also use a three-dimensional Fourier transform using a Tukey window in the 
vertical direction. The important observation here is that a power law scaling at small 
wave numbers is something that is characteristic of a small-scale turbulent dynamo, and 
is clearly different from the result of Kapyla et al. (20096) (see for example their figure 
9). If a large-scale magnetic field were to grow in addition to the small-scale one, this 
power law would not be observed. 

The lower two plots in figure [6] show the corresponding kinetic and magnetic energy 
spectra for cases R4a [Ta = 10^, A = 2) and R3a (Ta = 10'', A = 4). Although the aspect 
ratio is smaller in these cases than it is in the other two calculations that are illustrated 
in figure [6] the horizontal wave number has been defined in a consistent way across this 
set of spectra. This explains why the spectra that correspond to the smaller aspect ratio 
calculations cover a narrower range of horizontal wave numbers than those that are shown 
in the upper part of figure [6j Note also that the initial part of the kinematic dynamo 
phase has not been shown in these cases, which explains the absence of spikes in the 
magnetic energy spectra during the early stages of the calculations. Similar trends are 
observed in these smaller aspect ratio calculations. Again the kinetic energy spectrum 
peaks at smaller scales in the more rapidly-rotating case, whilst the magnetic energy 
spectrum always peaks at small scales. So, even in the most rapidly-rotating case, there 
is no evidence for large-scale dynamo action. It is more difficult to fit a power law to 
the low wave number range of the magnetic energy spectra in these smaller aspect ratio 
cases (which is why these have not been plotted), but the slope is not inconsistent with 
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what is found in the larger aspect ratio calculations. Similar results are obtained in the 
simulations with vertical magnetic field boundary conditions. 

These spectra strongly support the view that these calculations are all small-scale 
dynamos. This is consistent with the findings of previous studies of dynamo action in 
rotating Boussinesq convection ( Cattaneo fc Hughes|2006 Tobias et aZ.|2008j ) , so we can 
certainly rule out the possibility that the absence of a large-scale dynamo in these cases 



can be attributed to a lack of compressibility. Kapyla et al. (2009&) suggested that a 



necessary condition for the existence of a large-scale dynamo is that the Coriolis number 
Co ~ rid/Trf/inis exceeds some threshold value, which was found to be Co « 4 in their 
particular model. It is worth making the point that this definition of the Coriolis number 
may give a slightly misleading impression regarding the extent to which the convection 
is rotationally dominated, as it is defined in terms of the layer depth, rather than a 
characteristic length scale of the fiow (such as the horizontal integral scale, ^o)- As shown 
in the previous section, the integral scale is strongly dependent upon Ta, and is typically 
somewhat less than d. Hence it could be argued that the definition that is adopted by 



Kapyla et al. (2009&) actually overestimates the effective local Coriolis number. In any 
case, adopting their definition for the Coriolis number, our case R2b {Ta = 10^) has 
a value of Co w 0.4, R3c (Ta = 10^) has a value of Co w 6, whilst Co « 16 in our 
R4a (Ta = 10^) simulation, which is significantly higher than the approximate threshold 
value that was identified by Kapyla et al. ( 20096j . So although we cannot rule out the 
possibility that further increases in the rotation rate will make a significant difference, 
there is no evidence from these calculations to suggest that a large-scale dynamo would 
have been found by [Cattaneo fc Hughes^ (2006 J if they had simply increased the rotation 
rate in their simulations. 

As we have already noted, there are a number of differences between this model and 
the one considered by Kapyla et al. (20096 ), any of which may explain why they observe 
a large-scale dynamo. For example, they adopt a constant kinematic viscosity rather than 
a constant shear viscosity. Amongst other things, this implies that their Taylor number 
is independent of depth (it varies strongly with depth in the model that is described 
here). Given the perceived importance of rotation to the large-scale dynamo problem, 
we cannot rule out the possibility that this is the crucial difference. However, there is no 
hint of significant large-scale dynamo action in our quasi-Boussinesq case (R3d) , which 
suggests that this might not explain the differences between the two models. Another 
difference that is worth noting here is that the magnetic Prandtl number (the ratio of 
the magnetic Reynolds number to the kinetic Reynolds number) in the Kapyla et al. 



(20096 1 model was Pm w 2, which is rather less than the typical value of Pm « 5 for our 
simulations. As a result of this, particularly in their more rapidly-rotating calculations. 



Kapyla et al. (20096 1 were considering dynamos at comparatively low magnetic Reynolds 



numbers. It is therefore not inconceivable that ohmic dissipation might be playing a role 
in the formation of large-scale magnetic fields. Perhaps the most interesting possibility, 
however, is that the underlying stable layer might be playing a key role in the formation 



of large-scale magnetic fields. Indeed, Figures 4 and 7 of Kapyla et al. (20096) suggest 



that the strongest large-scale fields appear to accumulate preferentially around the lower 
boundary of the convectively- unstable part of the domain, which indicates that the stable 
layer may be an important part of their dynamo model. As mentioned in the Introduction, 



Tobias et al. (2008) found no evidence for large-scale dynamo action in their model of 



Boussinesq convection with an underlying stable layer. However, this may well be due 
to the fact that their layer was rotating less rapidly than that considered by |Kapyla| 
et al. (20096), who also failed to find a large-scale dynamo at modest rotation rates. 



Increasing the rotation rate in systems of this type tends to reduce the penetration 
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Run A Resolution 6 Ra Ta Hi C,o 
R2 4 256 X 256 x 120 3 5 x 10^ lO'^ 0.0044 0.3 
R3 4 256 X 256 x 120 3 1.5 x lO'^ lO'^ 0.0025 0.2 

Table 2. Summary of the parameter values for the simulations that were used for the 

measurement of the a-effect. 
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Figure 7. The mean electromotive force in the direction of the imposed magnetic field for (a) 
Ta = 10^ and (b) Ta = lO'^. The spatial average is carried out over the whole nume rical domain 



The value of Oimp is derived from the mean electromotive force using equation (5.141, and is 
normalised by ao = L'rms/3. 



depth of convective downflows overshooting into the stable layer ( Brummell et aL|2002 ), 
but why this should be beneficial for large-scale dynamo action remains an open question. 
As mentioned by Tobias et al. ( 2008 1 in their conclusion, the lack of magnetic buoyancy 
in their Boussinesq model might also be an important consideration in this context. 



5. Measuring the a effect 

As we have shown, despite there being a clear scale separation between the horizontal 
scale of the computational domain and the typical scale of convection, there is no evidence 
for significant large-scale magnetic field generation in our rapidly-rotating compressible 
convection simulations. Mean-field theory predicts that helical flows of this type should 
be producing an efiicient a-effect. If this indeed the case, the absence of any evidence for 
large-scale dynamo action is somewhat surprising. However, as described in the Introduc- 
tion, previous attempts to measure the a-effect in simulations of this type have yielded 
rather inconsistent results. In this section, we attempt to make direct measurements of 
the a-effect in our numerical simulations, focusing initially upon results that are obtained 
by using the standard imposed-field method. 
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5.1. Classical imposed-field method 

In the imposed-field method, the idea is to measure the mean (volume averaged) electro- 
motive force, £i, that is obtained when a mean magnetic field Bq is imposed across the 
flow. The corresponding value of a, denoted by aimp in the following, is then obtained 
from the following relation 

"imp = , (5.14) 

where Bq is the amplitude of the imposed mean field. We apply this method to the set 
of simulations that are defined by the parameters given in Table [2] Previous calculations 

suggest that an aspect ratio of A = 4 



(Hughes fc Cattaneo 2008 Kapyla et al. 2009a 



should be large enough to yield results that are not significantly constrained by the box 
size, so that is the aspect ratio that is adopted here. We have also carried out a few 
calculations with an aspect ratio of A = 8 in order to confirm that this is indeed the case. 
We consider two values of the mid-layer Taylor number, Ta = 10^ and Ta = 10^, as before 
choosing the Rayleigh numbers so that the global Reynolds number is approximately 150 
in each case. A uniform, steady horizontal magnetic field is imposed across the convective 
layer. This imposed field is aligned with the a;-axis, and we choose a field strength, Bo, 
that is very much less than the equipartition field B^q (typically Bq k, lO^^Bgq). This 
implies that the initial magnetic field does not exert any dynamical influence upon the 
flow. Given the difficulties that are associated with determining the value of a in the 
presence of a small-scale dynamo (see, e.g. Cattaneo fc Hugh"es||2009 1 , we choose values 
of Co that are large enough to ensure that these simulations do not exhibit small-scale 
dynamo action. 

As the magnetic fluctuations develop, we measure the a;-component of the mean elec- 
tromotive force £ = {ux B) across the layer. Figure [t] shows the evolution of aimp (as 
defined by equation (5.14)) for cases R2 and R3 (see Table[2]). Note that we average over 
the entire spatial domain in order to evaluate the mean electromotive force in our com- 
pressible model. This differs from the Boussinesq case that was considered by |Cattaneo| 
& Hughes (20061, where the mid-plane symmetry implies that only half of the domain 



can be used for this calculation. Note that it could be possible in principle to find a better 
way to compare Boussinesq and compressible simulations (by averaging about the zero 
in the relative helicity profile instead of averaging over the whole domain for example), 
but this approach led to the similar results in our particular case. Following [Cattaneo &:| 
Hughes ( 2006 1 , we also plot the cumulative average aimp defined by 



AT) 



T 



B~n 



dt 



(5.15) 



To normalise our results, we use the following isotropic expression derived from the first 
order smoothing approximation (see, for example, Kapyla et a/.|[2009a ) 

ao = ^C/rms ■ (5.16) 

For Ta = 10^, Uq = 0.1 whereas for Ta — 10^, ag — 0.05. In each case, figure [t] shows 
that the value of the mean electromotive force is strongly fluctuating with time, with a 
very small, poorly-defined time-averaged value. The local maxima of aimp are however 
of the same order as ao in both cases. These results are qualitatively similar to those 



obtained by Cattaneo & Hughes (2006) for their Boussinesq model. 

As shown in figure [2|Jc), the mean helicity changes sign within the domain. Therefore 
it is probable that some cancellation effects occur when the vertical average is taken in 
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Figure 8. The mean electromotive force in the direction of the imposed magnetic field for (a) 
Ta = 10^ and (b) Ta = lO''. The spatial average is carried out over the horizontal coordinates 
only. The thick line corresponds to the time average whereas the grey area is delimited by the 
peak values that are measured during the simulations. The thin dotted lines correspond to a 
selection of instantaneous profiles. 



order to calculate the mean electromotive force. This might imply that the value of the 
horizontally-averaged electromotive force, £{z) = {u x B)^, at some given depth could be 
substantially larger than its vertical average. To investigate this issue, we also measured 
the x-component of £{z). Figure [s] shows the depth-dependence of this horizontally- 
averaged quantity at different times for Ta — 10^ and Ta = 10^. The red straight 
line represents the time-average in each case. Although strong fluctuations are observed, 
even for this depth-dependent quantity the time-average is always very small. The only 
noticeable deviations from £x{z) = occur near the upper and lower boundaries, where 
it is worth noting that £x{z) and the mean helicity have opposite signs. 

5.2. Imposed- field method revisited 



As noted in the Introduction, Kapyla et al. (2010o) make the point that a complication 



arises when using the imposed- field method to measure a depth-dependent value of a. 
This is due to that fact that, even though the boundary conditions imply that the total 
horizontal magnetic flux is a conserved quantity, the horizontally-averaged magnetic field 
will, in general, be a function of depth and time. Figure[9]illustrates this effect. This figure 
shows some of the temporal variations of the horizontally-averaged magnetic field in the 
X and 2/-directions for the R3 simulation (R2 exhibits similar behaviour). Even though 
the vertical averages of < >z and < By >z are invariant (taking values of i?o and 
respectively), it is clear that these horizontally-averaged magnetic fields are strongly 



fluctuating quantities. So, even if the relation ( 5.14 1 still holds on average, since the initial 
flux is conserved, the local mean magnetic field is often very different from its initial 
value (possibly even zero, for some values of z and t). The horizontally-averaged mean 
electromotive force, at a given depth and time, should perhaps instead be interpreted as 
being proportional to {B^) {z) rather than Sq- To find a depth-dependent value for a 
at a given instant in time, it would therefore be necessary to normalise this quantity by 
a value of the mean field that may differ substantially from Bq. Where the mean- field 
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Mean magnetic field < B, >: Mean magnetic field < By >- 

Figure 9. Horizontally-averaged magnetic fields (a) {Bx)^ and (b) {By) ^ at different times. The 
thick dotted vertical line represents the initial value. These results correspond to simulation R3. 



is locally equal to zero, even this revised normalisation will not produce a meaningful 
result. 

In order to investigate this issue further, we carried out an additional simulation. This 
simulation is identical to R3 {Ta — 10^) in most respects. However, instead of solving 

for the horizontally-averaged magnetic field, we artificially constrain these mean magnetic 
fields to be constant everywhere in the layer. Despite this constraint, the associated 
magnetic fluctuations are still able to evolve according to equation (2.6), and no resetting 



process is used (as in Kapyla et al. 2010a). This calculation will allow us to test the 



hypothesis that the strongly fluctuating a-effect that was measured in case R3 was due 
to the presence of strong fluctuations in the horizontally-averaged magnetic field. In the 
test-field method, the magnetic fields are initially imposed (usually with some non-trivial 
spatial dependence), but these large-scale test fields are normally kept constant over time 
whilst the resulting magnetohydrodynamic fluctuations evolve. So in some sense, what 
we are doing here is equivalent to the test-field method, but with a spatially uniform 
mean field. 

The results from this new simulation are shown in figure [TO] where the measured value 
for a is denoted by ckiorn. In figure 10 a), the thick lines correspond to the results of the 
modified simulation whereas the thin lines are exactly the same results that are shown in 
figure ^h). The full spatial average of the mean electromotive force is now well-defined 
with a clear positive value. This value is roughly equal to half of oq. As shown in fig- 
ure |10[ b), the depth variation of the induced electromotive force shows significant positive 
values for a for < z < 0.6 and smaller negative ones for 0.6 < z < 1. This is a mirror 
image of the variation with depth of the kinetic helicity, which is illustrated in figure 
|2jc). Note also that the induced electromotive force is still fluctuating, as expected from 
its turbulent origin, but the fluctuations are less pronounced now that the variations in 
the mean fleld have been neglected. The results from this modified imposed-field method 
are (at least qualitatively) consistent with results from previous studies in which the 
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Figure 10. Results from the simulation for which the mean induction equation ( 5.17[ ) is not 
solved, (a) The mean electromotive force averaged over the whole numerical domain versus time. 
The thin lines correspond to the results already presented on figure ^h). (b) As in figure [s] this 
plot shows the horizontally-averaged emf as a function of depth. The red line corresponds to the 
temporal average whereas the black thin lines correspond to instantaneous profiles at different 
times. 



a-efFect was determined using the test-field method. Therefore, by artificially suppress- 
ing the time-dependence of the horizontally-averaged magnetic field, we have produced 
test-field-like results with an imposed-field calculation. 

Having established that measurements of the a-effect depend crucially upon whether 
or not fiuctuations in the horizontally-averaged magnetic field are included in the model, 
it is important to consider whether or not it is reasonable to neglect these fluctuations. It 
is difficult to see how the classical imposed-field method could be producing an erroneous 
value for a when a volume average is used to define the mean field (because the mean field 
is uniform in space and constant in time), so that certainly suggests that the modified 
method is producing an incorrect measurement in this case. As always, the situation is 
more complicated when horizontally-averaged quantities are considered. Equation (5.14), 
and the more general expansion involving the mean magnetic field derivatives (1.2), as- 
sumes that there is a linear relationship between the horizontally-averaged electromotive 
force {u X B) ^ and the horizontally-averaged magnetic field, {B) ^. To test this assump- 
tion, we again applied the standard imposed field method to the rapidly rotating sim- 
ulation R3 in Table [2j Once this system had evolved to a quasi-steady state, we found 
the depth and time-averaged linear Pearson correlation coefficient between | {B) ^ \ and 
I {u X B) ^ I to be 0.1 ± 0.03. There is therefore very little correlation between the am- 
plitude of the mean field and the amplitude of the resulting mean electromotive force. 
A similar result was obtained for the dynamo case R4a (in which there is only a very 
weak mean field in the absence of an imposed magnetic field), where the corresponding 
amplitude correlation coefficient was found to be 0.05 ± 0.01. This lack of correlation in 
amplitude explains (at least partially) the complications that arise when trying to use 
the imposed field method to compute a depth-dependent value for a. Clearly the level 
of correlation between the mean electromotive force and the mean magnetic field will be 
higher when the fluctuations in the mean-field are neglected (as is implicitly the case in 
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the test-field method). So although we cannot rule out the possibility that this measured 
lack of correlation could be model specific, this result strongly suggests that our modi- 
fied version of the imposed-field method has produced an artificially large value for a. 
Furthermore, given the relationship between this modified version of the imposed-field 
method and the test-field method, we would argue that the test-field method would also 
tend to over-estimate a in this case. 

In summary, we would therefore argue that there is no compelling evidence for a 
well-defined, dynamically significant a-effect in this particular system. So although the 
imposed-field method also has its problems, the most likely explanation for the absence 
of a large-scale dynamo in these simulations is that the strongly fluctuating mean elec- 
tromotive force produces a negligible a-effect. 



6. Conclusions and discussion 

In this paper, we have considered several numerical simulations of dynamo action 
in rotating compressible convection. Various parameters, including the rotation rate, the 
aspect ratio, the extent of the thermal stratification, the magnetic Reynolds number of the 
flow, as well as the magnetic boundary conditions, were varied. In all cases it was found 
that there was some value of the magnetic Reynolds number at which it was possible 
to sustain a hydromagnetic dynamo. Mean-field dynamo theory suggests that rapidly- 
rotating convection (which possesses a high level of helicity) should be able to generate 
large-scale magnetic fields on a dynamical timescale. However, all of these dynamos 
produced an intermittent distribution of small-scale magnetic fields. Given that similar 
results were obtained by Cattaneo & Hughes (2006), we can conclude that the absence 
of large-scale dynamo action in their model cannot be attributed to the fact that their 
Boussinesq calculations neglect the effects of compressibility. At first glance, however, our 
results would appear to be inconsistent with those of Kapyla et al. (2009& ), who suggested 
that large-scale dynamo action should be possible in such systems provided that the 
rotation rate is sufficiently rapid that the Coriolis number exceeds some threshold value. 
Although the Coriolis number in our rapidly-rotating cases does exceed this threshold 
value, it should be stressed that there are a number of important differences between 
this model and that considered by Kapyla et al. (2009&). Perhaps the most important 
difference is the presence of an underlying stable layer in the calculations of [Kapyla et al.\ 
(20096). Certainly the large-scale magnetic fields in their model seem to be organised 
around the interface between the base of the convective layer and the top of the stable 
region, which suggests that the stable layer is playing a crucial role in the large-scale 
dynamo. Further study is needed to confirm this suggestion. 

Motivated by mean-field theory, we also measured the a-effect in our simulations. As 
was the case in the Boussinesq model of Cattaneo & Hughes (20061, the imposed field 
method failed to produce a dynamically-significant value of a (regardless of the rotation 
rate). Even at a magnetic Reynolds number that is below the threshold for small-scale 
dynamo action, this quantity exhibited strong fluctuations, with a cumulative average 
that was very much less than the rms velocity of the flow. Prompted by a remark made 
by Kapyla et al. (20106 ), we also considered the effects of artificially suppressing the evo- 
lution of the horizontally-averaged magnetic field. Having made this modification to the 
imposed-field method, we found a much larger value of a. However, we would argue that 
this modified imposed-field method over-estimates a because, by fixing the mean-field, 
it artificially increases the correlation between the mean horizontal magnetic field and 
the mean electromotive force (the amplitudes of which appear to be almost completely 
uncorrelated when the mean-field is allowed to evolve in time and space). A similar ar- 
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gument could be applied to the test-field method. Therefore we would suggest that the 
most plausible explanation for absence of large-scale dynamo action in this system is 
that these flows produce a fluctuating mean electromotive force with an inefficient (and 
probably poorly defined) a-effect. 

The other previous rotating convection calculation that appears to exhibit large-scale 
dynamo action is a rapidly-rotating, near-onset Boussinesq simulation that was carried 
out by'Stellmach & Hansen (2004). Although we have not yet carried out any simulations 
at comparable parameter values, we did investigate dynamo action in near-onset (low 
Rayleigh number) rotating convection, varying the Taylor number up to Ta = 10^. As 
would be expected, the hydrodynamic convective flow is fairly organised, and the kinetic 
Reynolds number is relatively small. When a seed magnetic field is introduced at low 
^0, the flow is destabilised, leading to an increase in the kinetic energy as well as super- 
exponential growth in the magnetic energy. Even in this case, however, we did not flnd a 
large-scale dynamo, with the magnetic energy spectrum peaking at a similar wave number 
to the kinetic energy spectrum. Stellmach & Hansen ( 2004 1 were able to consider a more 
rotationally-dominated parameter regime, and it is probable that this could explain the 
occurrence of large-scale magnetic fields in their model. In future work, we intend to 
explore this parameter regime in larger computational domains in order to understand 
the transition between this near-onset solution ( Soward|1974 Stellmach fc Hansen|2004" I 
and the more turbulent regime. 
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